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ABSTRACT 



NOAA-7 AVHRR data from April 1982 and 1983 were used to perform a global 
scale analysis of aerosol particle characteristics. Several improvements were incorpo- 
rated into an AVHRR multichannel satellite data technique developed by Pfeil (1986). 
This included better cloud and sunglint discrimination, removal of Rayleigh radiance and 
accounting for ozone absorption. The charactersitics analyzed were optical depth and 
Aerosol Particle Size Index ( 5 , 3 ). S ^2 provides the slope of the aerosol particle size dis- 
tribution curve. Both of these parameters were evaluated during several naturally oc- 
curring events, foremost of which were the 1982 El Chicon eruption and the 1982-1983 
El Nino-Southern Oscillation event. The results provided evidence that a significant 
amount of aerosol particles over marine regions are from land-derived sources. How- 
ever, the results also provided evidence that some marine aerosol particles may be of 
biogenic origins. 
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. INTRODUCTION 



It has been popular to "blame" weather patterns and their changes on a wide range 
of events. Some of these events are ridiculous, such as blaming floods on man's first 
landing on the moon, while others merit serious attention, such as an El Nino>Southern 
Oscillation and it's eflect upon global weather patterns. One source that has received 
serious attention in the last two decades is the eflect that aerosol particles are having 
on the global climate. 

In the middle 1960's, McCormick and Ludwig (1967) hypothesized that a buildup 
of anthropogenic atmospheric aerosol particles could increase the albedo, thus causing 
a long term cooling of the earth. Charlson and Pilat (1969) described how aerosol par- 
ticles aflect the energy absorbed by the earth-atmosphere system. In succeeding years, 
mounting evidence has indicated that atmospheric aerosol particles play a crucial role in 
the earth's environment. 

These particles have the potential to directly influence the earth's heat budget 
through absorption, scattering and emission of radiation. Atmospheric particles also 
have a role in atmospheric chemistry. They provide a direct eflect through chemical 
interactions with other atmospheric constituents, such as hydrogen sulfide, and by de- 
creasing incident solar radiation that aflects minor constituents such as ozone. Thus, 
aerosol particles can have a profound effect upon global ecology. 

Coakley (1976) presented a model of the earth's radiation budget that inferred a 
need for spatial distributions of aerosol particles. Meanwhile, Weiss et al. (1976) used 
a model by Chylek and Coakley (1974) in which actual aerosol particle measurements 
were compared with clear air assumptions. This clearly demonstrated a need for aerosol 
particle information within radiation budget models. Short term temporal changes in 
temperature were linked to changes in aerosol particles during a Soviet investigation of 
longwave radiation in the tropical Atlantic ocean (Zaitseva, 1976). They observed a 0.3 
to 0.4®C/day radiative cooling with haze from North Africa. 

In remote marine regions, far removed from continental influences, aerosol particles 
had been thought to be composed primarily of sea salts. Woodcock (1952 and 1953) 
measured sea salts in the marine boundary' layer and proposed that they were responsible 
for most of the cloud condensation nuclei in marine regions, Meinert and Winchester 
(1977) found that the distribution of aerosol particles over the ocean was bimodal and 
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postulated that this bimodality was caused by large marine particles (sea salts, 
supermicron size) and continental particles (land originated, submicron size). The col- 
lective work of Carlson and Prospero on Saharan dust during the 1970's documents the 
invasiveness of continental particles over marine regions (Carlson and Prospero, 1972; 
Prospero and Nees, 1976; Carlson, 1978; Prospero, 1982). 

Patterson et al. (1980) measured a clearly bimodal distribution in remote regions 
during Global Atmospheric Measurements Experiment of Troposheric Aerosols and 
Gasses (GAMETAG) flights. He postulated that the smaller sized particles were marine 
in origin, not continental. Direct observations made by USXS LYNCH and NOAA ship 
DISCOVERY in 1983 and 1984 (Hoppel et al., 1987) revealed peaks at O.Ol^m and 0.08 
to 1.5^m. The samples taken during these cruises revealed that the aerosol particles were 
totally marine in origin and that the particles were too volatile to be predominantly sea 
salts. Hoppel et al. (1987) has demonstrated a method whereby both peaks could be 
produced by a stepwise process that starts with Dimethylsulfide (DMS) and ends with 
sulfur based cloud condensation nuclei (CCN) 

Hoppel et al. (1987) produced the bimodal distribution by converting DMS, through 
photooxidation into a product of low volatility, methyl sulfonic acid (MSA). MSA 
forms particles through homogeneous nucleation. Condensation of MSA and water 
produces a surface area size increase. Condensation starts when the MSA particles 
reach a radius of 0.04^m. This was considered by Hoppel et al. (1987) as a threshold 
radius for cloud activation. The conversion of D.MS to MSA cloud condensation nuclei 
is a rapid process, usually taking less than 12 hours. In a bimodal distribution, 
nucleation forms the smaller size distribution while condensation resulted in the large 
peak. Fig. 1, from Hoppel et al. (1987), shows how DMS initiated aerosol particles can 
increase in size over a period of time to form CCN. 

Much of the early evidence for oceanic sulfur sources tried to show a link between 
the ocean and atmosphere through fluxes of biogenic hydrogen sulfide and sea salt 
sulfates (Kellog et al., 1972). Lovelock et al. (1972) was the first to postulate DMS as 
the prime oceanic sulfur source. Since then, more evidenced from Charlson et al. (1987) 
gives credence to Lovelock's original hypothesis. DMS also has been found to be 
strongly linked to algal productivity in a variety of settings and in tropical regions in 
particular (Andrea, 1985). Regardless of the source, oceanic sulfur is now postulated to 
contribute 30% of the atmosphere's sulfur. 

If oceanic sulfur does contribute significantly to the atmospheric sulfur budget, then 
oceanic changes affecting that sulfur source would indirectly alter the earth's heat budget 
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Fig. I. Observed evolution of particle-size distribution of DIMS irradiation, from 
Hoppel et al. (1987). 



through the production or nonproduction of sulfur-based aerosol particles. Charlson 
et al. (1987) proposed a biological feedback mechanism for the earth's radiation budget, 
illustrated in Fig. 2, that uses sulfur-based aerosol particles as the buffer. 

Almost the entire body of knowledge of aerosol particles over the ocean has come 
from direct observations by ships or aircraft. While providing excellent data, the area 
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Fig. 2. Biological feedback in the earth's radiation budget, from Charlson et al, 
(1987). 
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covered is limited. Also, this method of data collection does not provide for long term 
continuous coverage of a specific location nor does it give broad coverage of an ocean 
basin. 

Presently, the only method of continuously monitoring aerosol particles on a global 
scale is by satellite remote sensing. While unable to detect the type of aerosol particle, 
satellites measurements can infer the relative size (large or small) of marine aerosol par- 
ticles in a given volume of air. When correlated with other occasional in situ data 
sources, such as cruise data, they can give a good measuremnent of the spatial and 
temporal distribution of atmospheric particles. Satellite measurements of aerosol parti- 
cles made on a global scale can provide a better understanding of the complex interre- 
lationship between the atmosphere and the ocean. 

The major purpose of this research is to develop a global scale perspective of aerosol 
particles and their role in the earths environment. In this thesis, multi-channel satellite 
data are used to measure global scale patterns of aerosol particle characteristics. Im- 
provements in sunglint and cloud discrimination are incorporated into the procedure and 
the results are compared with the similar work done by Pfeil (19S6) on the same data set. 
Pfeil (1986) used the satellite-measured red and near infrared radiance to infer the size 
relationship of aerosol particles. This was accomplished by forming a ratio of the 
albedos, as determined from those two wavelenghts. 

Additionally, an improvement in optical depth determination over the oceans is ex- 
plored using a variable scattering phase function. A result of Pfeil's work is that marine 
aerosols are not uniform in size or composition. Therefore, a single, constant scattering 
phase function is inadequate to accurately describe the spectral scattering by aerosol 
particles. 

Finally, the beginnings of an interannual comparison is performed through analyz- 
ing data sets for the same month in two years. The data sets are from before and during 
the 1982-1983 El Nino-Southern Oscillation (ENSO) event. This thesis will show that 
variations in aerosol particle characteristics are observable during an event of the mag- 
nitude of an ENSO. 
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II. THEORY 



Aerosol particles are remotely-sensed by satellites through their interactions with 
various wavelengths of upwelled radiation. This interaction can be by absorption, 
emission or scattering. For the purposes of this discussion, only red visible and near in- 
frared wavelengths will be used. The reason for this limitation will be discussed below. 

There are three primary sources of upwelled radiation in cloud-free areas, scattering 
by molecules in the atmosphere, scattering by larger particles called aerosols and the 
earth's surface. In the red and near infrared wavelengths, aerosol particles are the 
dominant scatterer over the ocean. Except for sunglint, reflectance by the ocean is ver\' 
small and scattering by molecules does not vary spatially. 

Satellites detect the upward scattered radiation from aerosols particles as upwelled 
radiance. Durkee et al. (1986) showed that there is a direct correlation between the up- 
welled radiance measured at the satellite and aerosol particle optical depth. Since optical 
depth is linked to the size and distribution of aerosol particles it should be possible to 
relate the satellite detected radiance to the size distribution of marine aerosols. While 
it has not been possible to develop a satellite technique for directly obtaining an aerosol 
particle size distribution, Bullfinch (1986) developd an approach for infering the size of 
aerosol particles. He used variations in radiance detected at red and near infrared 
wavelengths for an aerosol particle as an indicator of it's size. This process can best be 
explained by understanding the processes involved in the radiant energy transfer through 
the atmosphere. 

A. RADIATIVE TRANSFER IN THE MARINE ATMOSPHERE. 

Liou (1980) gives the following as the radiative transfer equation (refered to here- 
after as RTE) for solar radiation in a plane parallel atmosphere. 

/xdL(0; 0)/d(5 = L{S, Q,) - {(DqIAt:) f L(S, Q')dQ' 

J 4n (1) 

-{u>al4n)nFap{a, 

where 

L= diffuse radiance, 

S = optical depth, 

cjo= single scattering albedo. 
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^ = cos 6 {cos6 - satellite zenith angle). 
fiQ = cos 6q {6q = solar zenith angle), 

Q. = solid angle(0,0) (</) = azimuth angle), 
p{Q., = scattering phase function and 

tiFq = incoming solar radiative flux. 

The first term on the right hand side of the equation accounts for attenuation 
through either scattering or absorption. The second term describes the intensity added 
to the beam by multiple scattering. The third term accounts for single scattering 
intensification. For a more complete discussion, see Liou (1980). 

Over the ocean the RTF can be simplified due to the radiative properties of maritime 
air. Whenever the optical depth is small, the average path length between particles is 
relatively large. In general, optical depth over the ocean is less then 0.2 (Pfeil, 1986). 
If the aerosol particles are assumed to be randomly distributed and there is relatively 
large distance between them, single scattering dominates over multiple scattering. 
Therefore, it is reasonable to assume that the multiple scattering term of equation (1) is 
very’ small and can be disregarded. At these wavelengths the upward radiated intensity 
from the surface is also small (Sellers, 1965). Thus the RTF, reduces to 

Z.(0; n, 4>) = + //o)W©)(l - (2) 

For an optically thin atmosphere, <5(1/^ + Ijpo) < < 1. Therefore, the RTF can be 
further simplified to 

L(0; fi- 4>) ^ (cooFo/4/.)/.(0)5, (3) 

where 0 is the single scattering angle. 

B. OPTICAL DEPTH 

Optical depth provides a measure of the mass of absorbing and emitting material in 
an air column. It is defined as 




0 



which is the extinction coefficient of the atmosphere integrated from the surface to the 
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height of the satellite. This is the total extinction of the atmosphere. The extinction 
coefficient is given by 



where nr^ is the cross sectional area of the particle and n(r)dr is the size distribution of 
the particles at that cross sectional area. is an extinction efficiency factor that varies 
with the radius (r) and composition, through the index of refraction (m). 

These equations show the relationship between optical depth and the aerosol parti- 
cle. /?„, is a function of the aerosol particle size and distribution. Since optical depth is 
a function of /?„, , it follows that changes in the aerosol particle size or distribution will 
be reflected in optical depth. 

Attendent with these changes in the aerosol particles is a change in the radiative 
properties of the atmosphere. Equation (3) show’s that any change in the optical depth 
results in a change in the diffuse radiance that leaves the atmosphere. These changes 
are not only important in the measurement of atmospheric properties, but they also 
produce effects that are widely felt by other users of the electro-optical spectrum. 

In addition to detected radiance, optical depth is a function of single scattering 
albedo (wq), satellite zenith angle (^), solar radiation (Fq) and scattering phase function 
(p(©)). Based upon earlier assumptions, the major variables in determining optical 
depth are satellite zenith angle and phase function. Durkee (1984) showed that satellite 
zenith angle is the most critical factor effecting the radiance-optical depth relationship. 
Phase function w'as determined to play only a small role. 

Phase function only plays a small role if the aerosol particles in a given volume of 
air have a constant size distribution. This is a good assumption to make for small scale 
analyses. It is a poor assumption for large scale or global scale analysis, w'here large 
changes in the aerosol particle size distribution w^ould be expected. When the aerosol 
size distribution changes, the phase function must change also. Fig. 3 illustrates how 
phase function can change with relative humidity. On a global scale analysis, changes 
in ambient relative humidity must be expected. Therefore, any large scale analysis, uti- 
lizing a constant phase function produces a significant error. A more accurate repre- 
sentation, albeit more difficult, is to calculate optical depth with a phase function that 
varies with the dominant aerosol particle. 




( 5 ) 
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Fig. 3. Marine phase functions (after Shettle and Fenn, 1979) 

C. AEROSOL OPTICAL DEPTH 

The satellite detected radiance from equation (3) is the sum of aerosol particle radi- 
ance and path-added Rayleigh radiance (L = -I- Path-added Rayleigh radiance is 

radiant energy that is due to Rayleigh scattering in the atmosphere. Rayleigh scattering 
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occurs when the wavelength of the radiation is much greater than the size of the scat- 
tering material. At red and infrared wavelengths, the scattering material is primarily due 
to atmospheric molecules. The Rayleigh radiance must be accounted for in order to 
obtain the radiance due strictly to aerosol particles. By subtracting the Rayleigh radi- 
ance, equation (3) then becomes 

La = L(0; f,; icooFoMp{Q)S^^ (6) 

where is aerosol particle optical depth. 

It is possible to express the optical depth of the atmosphere as a function of meas- 
ured diffuse radiance by solving for 

<5a-'(4m/cooLoW0)^a- (7) 

This equation provides an understanding of the radiative properties of the atmosphere 
and, as will be seen, an insight into the size and distribution of atmospheric aerosol 
particles. 



D. RADIANCE AND AEROSOL SIZE DISTRIBUTION 

It has been determined that the aerosol particle dependent terms in the RTE are 
larger at red wavelengths than at near-infrared wavelengths for small aerosol particles 
(Durkee, 1984). Since the diffuse radiance is proportional to aerosol particle dependent 
terms, spectral variations in radiance result from variations in the aerosol particle de- 
pendent terms. The spectral variation can be quantified by a simple ratio between red 
wavelengths and near-infrared wavelengths. This results in 

(^A)red (<^oFo|4^‘)p(Q)iS/Or^i 

(LJ„„ ~ (a>o/-o/4M)p(0)(^A)mr ‘ ^ ’ 

Since for small aerosol particles, the aerosol particle dependent terms are larger, the ratio 
will be larger when smaller aerosol particles dominate in an area. 

Earlier it was noted that there was a linear relationship between aerosol particle 
optical depth and satellite detected radiance. This relationship is shown in Fig. 4. This 
relationship is linear because oJq is approximately one at red and near infrared wave- 
lengths and L(0) changes are small with wavelength variation (Durkee, 1984). Fq varies 
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as a function of wavelength, but this only results in a proportionality constant. The 
value of ^ is the satellite view zenith angle, which does not change with wavelength. The 
ratio, called the Aerosol Particle Size Index (Sij), can now be simplified to 



^ (^A)red 

"" (^A)nir ' (-^aU ’ 

where is optical depth due to aerosol particles. 

Satellite derived measurements cannot provide the actual aerosol particle size dis- 
tribution curve, but the ratio described above does give the slope of that distribution 
curve. The steeper, the slope the greater the influence of small particles on the distrib- 
ution and vice versa. 
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Fig. 4. Satellite detected radiance vs. aircraft measured optical depth, from Durkee 
(1984). 
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HI. PROCEDURES 



A. DATA 

The data were collected by NOAA-7, a sun-synchronous polar orbiting satellite with 
an average altitude of 833 km. It's equator crossing time is approximately 0230 local 
standard time northbound and 1430 local standard time southbound. Upwelled radiance 
was measured by the Advanced Ver>' High Resolution Radiometer (AVHRR) in five 
spectral bands or channels. Table 1 provides the bandwidths of these channels. Chan- 
nels 1. 2, 3 and 4 are used in this study. 



Table 1. NOAA 7 ADVANCED VERY HIGH RESOLUTION RADIOMETER 
CHANNEL BANDWIDTHS, AFTER JENSEN (1986). 



Channel 


Bandn idth 


Radiance 


1 


0.58-0. 68jum 


Red Visible 


2 


0.725-1. 10.um 


Red Visible Near Infrared 


3 


3. 55-3. 93am 


Mid Infrared 


4 


10.5-1 1.3/nn 


Thermal Infrared 


5 


1 1.5-12. 5^m 


Thermal Infrared 



The source of AVHRR data was a multi-year archive maintained at the National 
Center for Atmospheric Research (NCAR), Boulder, CO. From the AVHRR data for 
April 1982, two six-day periods, 5-10 April and 20-25 April, were studied. The area 
covered in this subset was the North and South Pacific Oceans, from 50°N to 50®S and 
from 110®E to 70°W. In addition, AVHRR data for April 1983 were also processed 
from the NCAR archives. The April 1983 data provided full global coverage. These 
data were for the full month, but only the odd numbered days of the month were studied. 
April 1982 data were analyzed in the two time periods, the 5th- 10th and the 20th-25th, 
to provide a direct comparison with the work done earlier by Pfiel (1986). 

B. DATA PROCESSING 

The April 1982 data were processed at the Interactive Environmental Digital Anal- 
ysis Laboratory (IDEA LAB) located at the Naval Postgraduate School (NTS), 
Monterey, CA. The heart of this lab is a VAX 8250 computer, with graphics display 
terminals. The April 1982 data, which were global data set, was processed on NCAR's 
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Cray XMP supercomputer with the assistance of the NCAR Scientific Computing Di- 
vision. The April 19S2 data were digitally transmited to the NFS IDEA LAB when 
processing was completed. 

The IDEA LAB was used to produce visual images from the processed data. The 
ability to enhance these images while working at the IDEA LAB's work stations facili- 
tated the interpretation and analysis of the images. 

C. DATA ANALYSIS 

The data analysis scheme was based on that used by Pfiel (1986), with several im- 
provements. Additionally, several other atmospheric features were analyzed, extending 
the early work. The final processing produced the following data sets: 

Channel 1 Channel 2 radiance ratios. 

Optical depth using a marine phase function. 

Optical depth using a rural phase function. 

Optical depth using a variable Heney-Greenstien phase function. 

Channel 1 reflectance for low clouds. 

Channel 3 reflectance for low clouds. 

Number of partly cloudy pixels. 

Number of high cloud pixels. 

All of the results were stored and averaged in 1° by 1® boxes. The number of AVHRR 
pixels stored in each box varied, depending upon a series of pixel tests. Although cloud 
data were calculated, it was not evaluated as part of this thesis. When the final images 
were produced, each data point is one of these 1° by 1® boxes. 

In Pfiel (1986), sunglint was removed by eliminating the entire western half of each 
image from 40®N to 10®S. This approach was overly conservative. While it ensured that 
there would not be any contamination due to sunglint, it also removed data points that 
were not contaminated. Since sunglint is a function of viewing geometrx', it is feasible 
to remove it more carefully by calculating an elliptical-shaped area. This is smaller than 
the rectangle used in the earlier work. This approach added about 25% of the removed 
area back into processing, but still provided assurance that little or no sunglint con- 
tamination occurred. 

Cloud removal was changed to a purely objective analysis vice a subjective/objective 
analysis used by Pfeil. All the image pixels w^ere subjected to the following tests in the 
order listed below; 
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1. Test for bad pixels; If channel 2 or channel 4 radiance values are less then zero, 
after Rayleigh radiance was removed, then the pixel is classifed as "bad" and not 
processed. 

2. Test for high clouds: A pixel was considered to contain high clouds if it's channel 
4 brightness temperature was less then 273®K. 

3. Test for low clouds: A pixel was considered to contain low clouds if it's channel 2 
albedo was greater then 40%. 

4. Preliminar}’ test for clear pixels: A channel 1 /channel 2 albedo ratio was calculated 
and if that ratio was greater then 1.5 the pixel was tentatively classified as clear. 

5. Final test for clear pixel: The channel 2 albedo of the pixel tentatively identified 
as clear was compared with the four pixels above, below, to the right and to the 
left. All the pixels must fall with in five counts of each other for the pixel to be 
classified as clear, otherv^’ise it was counted as partly cloudy. A count is equivalent 
to 0.1069% albedo. 

Equation (3) expresses a relationship for the measured radiance in both channels 1 
and 2 over the ocean. In Pfiel's earlier work, total radiance was assumed to be domi- 
nated by aerosol particle scattering, therefore Rayleigh radiance need not be accounted 
for. However, this is not the case for optical depth. 

Preliminarx' tests indicated a large difference in optical depth for cases with and 
without Rayleigh radiance removed. For cases that included Rayleigh radiance, optical 
depth estimates averaged 0.22 higher than in cases where Rayleigh radiance was re- 
moved. This is a significant value since the highest optical depth estimate obtained was 
0.S2. The cases including Rayleigh radiance showed optical depth ranges from 0.18 to 
0.82. Optical depth in the cases where Rayleigh radiance was removed ranged from 0.06 
to 0.54 with an average value of 0.24. For this reason, path-added Rayleigh radiance 
was removed. 

The removal of the Rayleigh radiance also produced another change in the proce- 
dure. Previously, the ratios were calculated using the albedo from each channel, rather 
than the radiance. This choice was made based on computational considerations. To 
remove the Rayleigh radiance also requires the determination of the ratios using radi- 
ances. 

A total of three optical depth calculations were performed using three separate 
phase functions. As discussed earlier, due to the changing nature of aerosol particles 
over various oceanic regions, any one invariant scattering phase function was felt to be 
inadequate. The three phase functions used were a "classical" marine phase function 
(Shettle and Fenn, 1979), a continental rural type phase function modified for over 
ocean usage and a varying phase function. 
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The two constant phase functions were used to determine whether accurate results 
were obtainable using a simple phase function. Optical depths are much easier to derive 
using a constant phase function. However, the results may not be adequate for a global 
scale analysis. Using all three phase functions would provide a method of choosing the 
most appropriate one for global scale estimates of optical depth. 

Based upon the considerations just discussed, an appropriate phase function was 
chosen. Fig. 5 shows a composite of optical depths taken from several successive sat- 
ellite passes in a single day (5 April 1982). From this figure it can be clearly seen that 
the optical depth estimate that included the constant marine phase function (Shettle and 
Fenn, 1979) was inconsistent. The optical depths varied across the width of the satellite 
pass, from west to east, with a maximum in optical depth towards the west side, de- 
creasing towards the east. In contrast, optical depths are more uniform when a var\'ing 
single scattering phase function is used. The optical depths are fairly uniform across the 
width of the passes and the overall optical depth results were in the expected range of 
less then 0.25. The optical depths from the invariant phase function were much higher 
than expected, with large areas having values above 0.5. 

A problem with using an invariant phase function is that it is based upon a volume 
of air containing a uniform aerosol particle size distribution. Any significant deviation 
from that aerosol particle size range and the scattering angles predicted by that phase 
function will be erroneous. A major premise of this work, based on the results of Pfeil 
(1986), is that the aerosol particles over the ocean are not uniform in size or distribution. 
It is obvious , from the small sample shown in Fig. 5, that the constant phase function 
is resulting in optical depths that are too high, greater than 0.5, (ie., greater than 0.5), 
on the western side of the images. 

The relationship between optical depth and phase function was presented in 
equation (7) from Chapter II. Optical depth, is inversely proportional to the p(0). 
Thus, when the western portion of the satellite images show evidence of optical depths 
too high, p{Q) is too small. Since the overall optical depths are too high , the constant 
phase function is most likely predicting too much forward scattering for the aerosol 
particles present. From Mie scattering theor>% it can be deduced that the aerosol parti- 
cles actually present are smaller then the aerosol particles from which the marine phase 
function was developed. 

The variable single scattering phase function is superior because it starts to account 
for differences in the aerosol size distribution. It works on the basis of choosing the best 
phase function, based upon an inferred aerosol particle size. By providing the correct 
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Fig. 5. A comparison of optical depths as determined from an invariant and a vari- 
able phase function. (5 April 1982) 
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phase function for an aerosol particle distribution, a more accurate optical depth esti- 
mate is produced. 

The relative sizes of the aerosol particles are determined by the 5,2 value and based 
upon this ratio a phase function is calculated. This provides for a close fit between 
phase function and aerosol particle size. In Fig. 6, the higher ratios, which are due to 
smaller aerosol particles, have more backscattering radiance than the lower ratios, which 
are due to larger aerosol particles. This provides the correct P{Q) needed to estimate the 
optical depths expected from small particles. 

The variable phase function used was a two term Heney-Greenstein scattering phase 
function (Lenoble, 1985). 



7 ^( 0 ) = 



(i+ 5 ?- 2 g, cos(e)) 



(I - »0(i-g|) 

(l+?2+2g2COS(0)) ’ 



(10) 



where IF is a weighting factor and and g 2 are assymetrv' curve factors. All three of 
these factors are functions of 5,2: 



IF= 1.2 -(0.25x5,2), 


(11) 


= 1.04 — (0.08333 X 5j2)a 


(12) 


^2 = 1.2 -(0.58333 x 5i2). 


(13) 



As equations 10-13 show, the phase function varies as a function of 5,2 . The basic 
assumptions are that in areas that have relatively low ratios, the aerosol particles are 
large and in areas of high ratios the aerosol particles are small. Larger aerosol particles 
are assumed to follow a marine type phase function, while the smaller particles are as- 
sumed to follow a rural type phase function. As the ratios change, the factors used 
change, thus varying the phase function from one similar to a marine phase function to 
one similar to a rural phase function. Fig. 6 shows some of the variable phase functions 
calculated along with the 5,2 used in the calculations. 

One final change in the the procedure is an accounting of the effects of ozone. 
Ozone is a weak absorptive substance in the red and near-infrared wavelengths. It's ab- 
sorption is greater in the red visible spectrum than the near infrared. The major ab- 
sorption is at a wavelength of 0.69^m (Fleagle and Businger, 1980), which is found in 
channel 1 of the AVHRR. A much smaller amount of radiative absorption occurs in 
channel 2. This produces an unequal amount of ozone absorption of both incoming and 
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Fig. 6. Aerosol Particle Size Index and associated phase functions. 



outgoing radiation. Since the amount of loss due to ozone is a function of wavelength 
and view zenith angle, it must be properly accounted. Although ozone absorption is the 
smallest of the changes, it can produce up to a 10% change in estimate of optical depth. 
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IV. RESULTS AND DISCUSSION 



A. AEROSOL PARTICLE SIZE INDEX (CH1/CH2 RATIO) AND OPTICAL 
DEPTH 

I. Intercomparison of Chl/Ch2 ratio techniques 

First, an analysis was performed comparing the current data processing scheme 
with the results obtained in Pfeil (1986). A direct comparison is only possible in the ratio 
results, because Pfeil (1986) calculated albedo and did not estimate optical depth. 

Fig. 7 and show the S ^2 values from April 1982 for Pfeil's earlier work and the 
current study. Both results agree in a maximum Sjj value in the tropics, with the highest 
values located close or over the equator. They both indicate that the S ^2 values decrease 
poleward of the tropics, but with a big difference between them in the amount of de- 
crease in the northern hemisphere. Pfeil's results indicate that the decrease in the north- 
ern hemisphere is equivalent to the decrease in the southern hemisphere. His calculated 
Si 2 values ranged from a maximum 2.0 at the equator, decreasing poleward to an average 
of 1.2 in the midlatitudes. Effectively, the northern hemisphere is equal to the southern 
hemisphere. The present work does not have equivalent decreases. The average S ^2 in 
the northern hemisphere for the present study is 1.9, while in the southern hemisphere 
it is 1.65. 

This comparison also indicates that Pfiel's 1986 results underestimate the S ^2 
globally. His maximum in the tropics averages 1.8, which is lower than the present 
work's tropical average of 2.2. Pfeil's S ^2 estimates also are too low in the mid-latitudes, 
as discussed above. 

On a smaller scale, there are dilTerences in the Gulf of Mexico, the northwestern 
coast of Australia and the southern coast of Australia. Pfiel's study, shows a minimum 
in 5 j 2 values in the Gulf of Mexico. In the present work, the values observed in the Gulf 
of Mexico are much higher, approximately double Pfeil's results. Off the northwestern 
coast of Australia, Pfiel's work showed some of the highest values observed, greater than 
2.0. Also, off the southern coast of Australia is a relative maximum, approximately 1.8. 
While the S ,2 values observed in both studies for these areas were approximately the 
same, they are not maximums in the present study. 
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There are several reasons for these differences. These include the removal of 
Rayleigh radiance and accounting for ozone absorption in the present scheme along with 
a better cloud discrimination technique. 

Rayleigh radiance removal has a two-fold influence. The first is that it lowers 
the channel one radiance more than the channel two radiance. This produces a final 
Si 2 value that is lower than the produced with Rayleigh radiance included. Also, the 
removal of Rayleigh radiance eliminates most of a geometric dependence found in the 
earlier 5i2 values. In Fig. 7 a strong dependence on geometry’ is visible, particularly in 
the late April image. Higher values are observed towards the eastern half of the passes. 
This is due to the satellite sun-geometr\’, where the incoming incident solar radiation, 
because NOAA-7 is an afternoon satellite, is along the western edge of the pass. This 
produces a maximum in path-added Rayleigh radiance along the eastern portion of the 
satellite pass, along with a greater amount in channel 1 than in channel 2. This most 
likely accounts for some of the higher S ^2 values observed in the northern hemisphere, 
since the solar elevation is higher in the northern hemisphere at this time. 

Another contributing factor to the northern hemisphere having lower Sj, values 
in the earlier study is that the cloud discrimination technique is not as effective. Clouds, 
with their large water droplets, would produce low S ^2 values. The blank area in 
Fig. 8. which is an area that does not contain any clear pixels, has almost a sixfold in- 
crease over the blank area in Fig. 7. This extra blank area can be equated to an increase 
in the ability of the current method to eliminate cloud contaminated pixels. The earlier 
approach allows a number of cloud contaminated pixels to affect the albedo ratio values. 

When the maximum S ^2 values are compared between the two sets of results, it 
should be noted that the current work produces higher maximum values than Pfeil's 
study. This is probably due to ozone absorption correction. This absorption varies with 
wavelength and path angle through the stratospheric ozone layer. The longer near- 
infrared wavelength is absorbed less then the shorter wavelengths. The result is some 
elevation of the S ^2 values with an ozone absorption correction. 

The differing sunglint removal techniques has two effects. The first effect is that 
the artificial line that is observable in the tropics in Fig. 7 is no longer visible in the later 
work (Fig. 8). This line, that was produced previously, was strictly due to a straight 
east-west line being the northern and southern most limits of the sun glint removal re- 
gion. The second effect of the sunglint removal technique is to add more useful data 
pixels into the statistics. This likely accounts for the smoother pattern shown in 
Fig. 8. 
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Fig. 7. Aerosol Particle Size Index for early and late April 1982 from Pfeil (1986). 







Fig. S. Aerosol Particle Size Index for early and late April 1982. 
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2. Albedo and Optical Depth 

Pfeil (1986) did not estimate optical depth. Instead albedo was used to estimate 
atmospheric conditions. It is much better to use optical depth, which gives an under- 
standing of the transmittance capabilities of the entire atmospheric layer, rather than 
albedo which measures the total amount of radiance reflected back to the satellite. From 
Fig. 9, which reproduces Pfeil's albedo results, it is easy to see that albedo is affected 
by the sun-satellite view geometry. This may be due to the effects of Rayleigh radiance, 
which was discussed above in connection with albedo ratios. 

Fig. 10 shows the optical depths obtained using the present method. In this 
figure, a high optical depth region extends west into the North Pacific Ocean from 
Central America. There is a similar feature in Fig. 9. However, in Fig. 10 the feature 
is a continuous discrete band that has the appearence of being a natural occurrence, 
while the feature observed in Fig. 9 is very disjointed and looks as if it may be an artifact 
of the processing. The optical depth values, although high, are consistently in the 0.3 
to 0.4 range, without many discrete breaks. The albedo starts at 12% and decreases to 
2%, then jumps back to 12% as the feature is tracked westward. There is a similarly 
occurring feature in the later half of April 1982 in the western North Pacific that shows 
the same differences as those just described for the eastern North Pacific. 

Cloud contamination appears to be more evident in albedo than optical depth. 
In the albedo results, high values were located in all areas normally associated with large 
cloud coverage. In particular, cloud contamination is the cause of high albedos in the 
region of the South Pacific Convergence Zone. In the same region, the optical depths 
were all low. While this area is a cloudy region, the cloud discrimination removed all the 
cloudy pixels when computing optical depth, while the albedo technique did not perform 
as well. 

3. Aerosol Particle Size Index and Optical Depth 

Analyses were performed intercomparing early April 1982, late April 1982 and 
April 1983. As discussed in Chapter II, these time periods were not identical. The early 
and late April 1982 data comprised two separate six-day periods, while the April 1983 
data consisted of the 15 odd numbered days, covering the full month. Additionally, the 
1982 data were primarily in the North and South Pacific and the 1983 data was global. 
However, these differences did not pose any difficult problems in the analysis. 

In April 1982 the dominant synoptic pattern was the slow development of a 
strong long wave ridge in the mid North Pacific Ocean, with a long wave trough axis 
over eastern China. Early in April 1982 the polar jet south of Japan was primarily 
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Fig. 9. Albedo for April 1982 from Pfeil (1986). 
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Fig. 10. Optical depth from April 1982. 
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westerly. With the development of the long wave pattern described above, the jet turned 
towards the north slightly, as southwesterly How. At the same time, strong cyclonic 
activity, associated with the long wave trough, over the Gobi desert was causing strong 
dust storms. These conditions resulted in favorable conditions for long range dust par- 
ticle transport over the North Pacific. 

The major aerosol particle producing event, for 1982, was the eruption of the 
El Chicon volcano (17.3“N, 95. 2® W). There were several eruptions between 28 March 
and 6 April 1982. Debris from these eruptions moved both west and east (Krueger, 
1983). This would produce some significant effects upon the optical depths and the S ^2 
values. 

A basic overall pattern in optical depth estimates is that there are lower values 
in the southern hemisphere compared to the northern hemisphere. A typical value in the 
southern hemisphere is 0.1 or less, while in the northern hemisphere, optical depth is 
consistently 0.15 or greater. This is shown in Fig. 10, April 1982 optical depth estimates, 
and in Fig. 11, which presents the April 1983 optical depth estimates. This optical depth 
pattern is a result of the greater land mass located in the northern hemisphere. Land- 
derived aerosol particles are in a much larger quantity in the northern hemisphere, both 
naturally occurring and man-made, which results in the larger optical depth estimates in 
the northern hemisphere. 

In early April 1982 (Fig. 10) there is a maxima greater than 0.25, found off the 
west coast of Mexico that and extending west about one third of the way across the 
North Pacific Ocean. Other maxima of similar values are located along most of the 
eastern Asian coast line in early April, along with a few isolated areas in the mid North 
Pacific. In the Gulf of Mexico there is a realtively high maxima (almost 0.5) that extends 
north and east. 

In the later half of April 1982 many of the same maxima exist, with one major 
exception. The maxima off the coast of China in the East China Sea extends far to the 
east and has increased in thickness in some locations to 0.6 or greater. It appears to be 
merging into the maxima from the west that is still evident. 

The high optical depths just noted are, with a few exceptions, all in areas which 
are adjacent to large land masses. Two of these maxima, along the East China Sea and 
the southern Mexican coast, are located near major aerosol particle producing events. 
Gobi desert dust and El Chicon debris would be expected to produce high optical depths. 

April 1983 shows evidence of dramatic changes in optical depth from the previ- 
ous year. The strong maxima that extended across the North Pacific Ocean is no longer 
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Fig. II. Optical depth from April 1983. 
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present, although there is a very small area that has very high values of optical depth, 
greater than 0.25, along the western coast of Central America. There are also high val- 
ues observable, again, in the East China Sea and there is a region of high optical depths 
surrounding the island of Borneo. These results are consistent with the 1982 results. In 
1983, El Chicon was not erupting, so optical depths would be expected to be lower in 
that region. The maximum near Borneo coincides with one of the worst series of bush 
fires ever recorded in that part of the world (Voice and Gauntlett, 1984). 

Elsewhere around the globe in the 1983, there is a ver>' high maxima extending 
out of western Africa into the North Atlantic, with thickness values exceeding 0.5 close 
to the coast. There are also maxima greater than 0.25 in the Amazon delta basin, the 
Gulf Coast of the United States and near dense coastal population centers. The optical 
depths also appear to be reasonable over smaller bodies of water, which was not ex- 
pected. 

In April 1982 the dominant Sj 2 feature is high values, above 2.5, in the tropical 
North Pacific Ocean (Fig. 8). The values were higher in the eastern ocean averaging 2.4, 
with a decrease westward to an average of 2.2. As noted above, the S ^2 values decreased 
poleward, with a greater decrease in the southern hemisphere vice the northern hemi- 
sphere. 

Also, in early April 1982, there is a maximum 5,2, greater than 2.5, along the 
northern coast of Australia which is not observed later in the month. There is also a 
secondary' diffuse S ^2 maximum, aproximately 2.25, spread over a large portion of the 
Philippine Sea in the early part of the month. This changes dramatically in the second 
half of the month, when a strong 5i2, ranging from 2.25 to 2.50, comes off the southern 
coast of China, extends through the East China Sea and east into the Philippine Sea 
before curving to the northeast in to the mid north Pacific Ocean. Other 5,2 maxima 
that are observed in both early and late April 1982 are along densely populated coastal 
areas, such as the Chinese coastline, around Japan and the west coast of the United 
States. 

Just barely visible along the eastern edge of Fig. 8 is the Gulf of Mexico and the 
Caribbean. Careful note should be made of the ver>’ high 5,2 located to the south of 
Florida, as well as the one found along the Gulf Coast of the United States. 

April 1983 (Fig. 12) shows a dramatic change in the overall 5,2 pattern observed 
in 1982. The tropical maximum that extends across the North Pacific is dramatically 
diminished. Where 5,2 values were 2.0 or greater in 1982, they almost never exceeded 
2.0 in 1983, and were often much less than that. The strong maxima found along the 
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west coast of northern South America is just barely visible and is much weaker. Along 
the coast of the western United States and Mexico the values are greatly diminished with 
the exception of isolated maxima near Los Angeles and San Francisco, which retain S 12 
values greater than 2.5. 

The pattern of decreasing S 12 values poleward, with a greater decrease in the 
southern hemisphere, continues in April 1983. Since 1983 is a global scale analysis, this 
pattern is observed throughout the entire globe as can be seen in Fig. 12 

The maxima south of Florida noted in the 1982 analysis appears again in 1983. 
Again, it is one of the highest S ^2 maxima found around the globe, maintaining a value 
in excess of 3.0. Maxima along the Gulf Coast of the United States are still visible as 
well as coastal maxima along the Chinese and Japanese coastlines. Elsewhere around 
the globe, maxima are found in coastal regions near dense population centers, such as 
Bombay and Calcuta as well as in some areas of high industrial output, like the oil re- 
fmer\’ areas of the Persian Gulf One interesting region of high 5^2 values is the broad 
delta region of the Amazon river. 

Due to the nature of requiring clear air in a low optical depth region, it was 
expected that this procedure would only work over oceanic regions. However, maxima 
are clearly visible in the Great Lakes region near Detroit, Chicago and Clevland. Results 
are also obtainable in other small bodies of water, such as the Mediterranean, Black and 
Caspian seas. 

The S 12 results can. in many cases be correlated to the optical depth results. 
High S 12 values were often located in proximity to land masses and naturally occurring 
events. However, some high values were observed in regions far removed from land 
sources, without any apparent means of reaching those remote areas. These will be 
discussed in more detail below. 

At this point it should be remembered that the S ^2 value is the slope of the aer- 
osol size distribution curve. It does not provide any absolute value as to the number 
of various sized aerosol particles present. The slope may remain the same for a large 
number of aerosol particles as a small number of aerosol particles. What it does reveal 
is the relative size of the dominant aerosol particle. High values denote small aerosol 
particles (submicron) and low values denote large aerosol particles (supermicron). 
Therefore high S 12 values indicate that small aerosol particles dominate, but there is no 
evidence for the number of aerosol pariticles present. 

In several areas, high S 12 values are collocated with high optical depths. This is 
seen in the eastern North Pacific along the coasts of Mexico and the United States as 
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Fig. 12. Aerosol Particle Size Index, April 1983. 
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well as along the eastern Asian coast line. This would indicate that there is a ver>‘ high 
number of aerosol particles in the atmosphere, hence the thick optical depth values. 
Since the transmittance of the atmosphere decreases as optical depth increases, high 
optical depths are indicative of some scaterring or absorbtive agents. The high S 12 values 
in the same regions show that the distribution of the aerosol particles are weighted to- 
wards smaller particles. Therefore, it is reasonable to conclude that large amounts of 
aerosols, continental in origin, are blowing ofT the land and out over water in these areas. 
They do not extend great distances, because they tend to settle out before they have a 
chance to travel extensively. The particles that settle out are also large enough to form 
CCN. As such, they will also settle from rainout. 

The maximum in the ratio extending out into the North Pacific Ocean from the 
El Chicon region is also collocated with a ver\- high optical depth value. This extension 
into the North Pacific is from strong low and mid level easterlies carrying the volcanic 
debris out over the ocean. Note however, that the extent of the high optical depths is 
not as great as the Sjj values. This is most likely due to the settleing of the low level 
debris out of the atmosphere as it moves away from the volcano. 

Krueger (1983) showed that the debris from El Chicon moved east in addition 
to moving west. He proposed that this eastward movement was from debris that had 
been ejected up into the stratosphere. High optical depths are observed extending north 
and then to the east of El Chicon in early April 1982. They are not evident in late April 
to the extent that they are in the earlier portion of the month. An analysis of the the 
upper-level winds. Fig. 13, shows southerly flow over the El Chicon region, which then 
turns into westerly flow further north. This appears to support the contention of 
Krueger (1983) that some of the debris blew to the east. 

In the western North Pacific, the high optical depths are in the same regions 
along the coastlines as the high S 12 values in the early part of April 1982. Then, in the 
second part of the month they both extend out into the North Pacific. Earlier, in the 
middle of the month, several strong cyclones swept through the Gobi desert. It is not 
uncommon for the dust from the Gobi desert to blow long distances out over the ocean 
(Shaw, 1980). The high 5,2 values in late April follow the mid- and upper-tropospheric 
wind patterns given by Fig. 14. Therefore these high 5,2 values are probably from dust 
originating in the Gobi desert. Notice, again, how the optical depths decrease with the 
distance from land. This also provides evidence that the source of the aerosol particles 
producing the high 5,2 values are settling out of the atmosphere as they travel. 
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Evidence of blowoff from Asia, most likely Gobi desert dust, was also seen in 
April 1983. High S ^2 values were all along the coast, with decreasing values further east 
out over the ocean. Higher optical depths were also noted along with these higher S 12 
values. The extent of this blowolT was not as great or as strong as the 1982 event. This 
may be due to weaker winds over the interior portion of the Asian continent in April 
1983 versus April 1982. 

Pfeil (1986) did an extensive analysis of the high Sjj values found along the 
northern coasts of Australia. While his results covered a much larger area, particularly 
off the northwestern coast Australia in the latter half of April 1982, they were still in a 
similar range as the present study. He proposed that strong offshore flow accounted for 
the increased 5i2 values over water. However, the optical depths, which he did not 
evaluate, show very’ low thickness values in this region, with a small exception in the 
Gulf of Carpentaria and south of Port Moresby. This means that there may be aerosols 
blowing off the coast of Australia, however, they must all be small and in a much lower 
concentration then seen in the Gobi desert dust or the El Chicon debris. Therefore, 
north of Australia, small aerosol particles dominate, but there are less aerosol particles 
overall than in the two cases discussed above. 

The tropical region extending across the North Pacific is similar to the areas 
north of Australia, in that there are high S ^2 values but low optical depths. This area 
had very high S ^2 values in the east, that slowly decreased westward. There is no similar 
trend with regards to the optical depths. This denotes that there are very high levels of 
small particles compared to large particles. However, the overall number of particles 
must be reduced compared to the Gobi dust and El Chicon debris. Pfiel (1986) had this 
same feature in his results and hypothesized that the convergence of equatorial easterlies 
might account for this maximum, with a continental derivation of aerosol particles. 
However, up to this point, the cases that have been obser\’ed containing continentally 
derived aerosol particles long distances from the coast start with very' thick optical 
depths along the coastline. This is not evidenced in this case. Also, the line of maximum 
values is to the south of the April ITCZ. There must be some other source, most likely 
an oceanic source. 

Andrea (1985) found that dimethylsulfide (DMS) concentrations are veiy high 
in upwelled tropical waters. Clarke (1987) found an increase in atmospheric sulfur 
compounds in the equatorial regions. In particular he found veiy^ high elevations of 
small sulfur based aerosol particles all along the equatorial upwelling regions. Given 
that coccolithophorids are prime producers of DMS (Andrea and Raemdonck, 1983) and 
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that the dominant phytoplankton in equatorial regions are coccolithophorids (Hasle, 
1959; Okada and Honjo, 1973). then it is conceivable that the high Sjj values in the 
equatorial regions are from a flux of biogenic sulfur across the air sea interface with a 
subsequent conversion to sulphate particles as hypothesised by Charlson et al. (1987). 

The first observation concerning the overall pattern in April 1983 is the lower 
values for the aerosol particle size index in the tropics. If as surmised above, that the 
equatorial maxima are due to marine biogenic sulfur compounds, then this reduction 
might be expected. The time period, April 1983, was during the 1982-1983 El Xino- 
Southern Oscillation event. Starting in October 1982, upwelled ecosystems in the east- 
ern equatorial Pacific began a transition from normally high productivity to low 
productivity that did not return until July 1983 (Barber and Chavez, 1983). Feldman et 
al. (1984) found a complete change from the normal phytoplankton distribution by 
March 1983 in the eastern equatorial Pacific Ocean. This tends to support the hypoth- 
esis of Charlson et al. (1987) that phytoplankton are responsible for a large portion of 
marine aerosol particles. 

Outside of the North and South Pacific Oceans, the most noticeable feature in 
April 1983 was the high optical depths observed extending from the w^est African coast 
into the North and South Atlantic. This is consistent wdth observations confirming 
Saharan dust from the coast of Africa extending into the western Atlantic basin and into 
the Caribbean (Prospero, 1982). 

The other interesting feature is the high Sjj values and thick optical depths 
found along the Amazon river delta. This region is known to be one of the most highly 
productive biological regions in the w’orld. Lovelock et al (1972) discussed the pro- 
duction of aerosols from terrestrial forests as a biological emission. This may possibly 
account for some of this increase. A second possibilty is the clearing of land in the 
tropical rain forest through large scale burning. This would also introduce large 
amounts of aerosol particles into the atmosphere. This is most hkely the cause of the 
high values observed on either coast of Central America in April 1983. 

The most consistent area of high S ^2 value is the location south of Florida that 
has these values both in 1982 and 1983. The consistency in position and strength is 
extrordinar\^ While there are several tropical islands that may provide a contamination 
from land into the process, this is not obserx'ed in other areas, such as the Haw'aiian Is- 
lands. In fact the process for determining 5,2 values should produce low'er values over 
land, which is evidenced in Fig. 8 and Fig. 12. These islands do not contain dense 
population levels nor do they have the heavw industry they may account for these values. 
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Neither are the waters surrounding the islands an^■^vhere near as productive as the 
equatorial upwelled regions. Thus this region remains an anomalous region that requires 
more study. 

B. CLEAR PIXEL POINTS 

The number of clear pixel points was recorded in order to determine the statistical 
relevance of the data set. Additionally, the converse of the clear pixel points would 
provide a indication of the cloud representation during the month. In Fig. 15 and 
Fig. 16, which show the pixel counts for April 1982 and 1983, respectively, there is a 
maximum in clear pixel points in the tropics, with 1983 shifted slightly northward of the 
1982 position. 

As can be seen from Fig. 15 and Fig. 16 there is consistently more than 500 pixel 
points in most of the V by V boxes in the images. This means that the averaging 
technique used over 500 data points for each S ^2 value and optical depth estimate. In 
many cases the number of data points exceeded 1,000 points. 

One thing that can be noticed is that the number of points appears to increase in the 
center of the passes. This is due to the overlapping nature of the passes and the removal 
of pixel points through sun glint removal or the elimination of pixels from the edges 
because of reduced resolution (16 pixels from each edge of a pass). 

In the tropical area of very high Sj 2 values during April 1982, there is a corre- 
sponding high number of clear pixel points. This would seem to indicate that these ver\’ 
high Sj 2 values, with low optical depths, do not produce a significant amount of clouds. 
Just to the north of this high S ^2 area is the ITCZ. Notice that the S ^2 values decrease 
as well as the number of clear pixels. The easterly winds that converge at the ITCZ can 
move these small aerosols, which will be undergoing nucleation and hydroscopic growth 
with time, into that ITZC. By the time they reach the ITCZ, they will be large enough 
to form cloud condensation nuclei. Hoppel et al. (1987) demonstrated how DMS, the 
sulfur compound that may be the cause of these aerosol particles, can undergo the 
neccessary chemical and physical growth processes in order to reach cloud condensation 
nuclei size. This also provides more evidence that the tropical aerosol particles may be 
biogenic. 
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Fig. 15. Number of dear pixels used for April 1982. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Several conclusions can be drawn from the results discussed above. These can be 
split into two major sections. The first section deals with the technique and it's im- 
provements. The second focuses on the results of a satellite-derived global analysis of 
optical depth and Aerosol Part Particle Size Index (Sjj). 

It is obvious from the results that the improved techniques produced some signif- 
icant results. Foremost among these was the introduction of a variable phase function. 
Without the variable phase function, it would not have been possible to perform a global 
scale averaging of optical depth. All the previous work estimating optical depth indi- 
cated that, while var\’ing the phase function was desirable, it was considered too difiicult 
to implement. This work showed that not only was it essential to var>’ the phase func- 
tion, but that it could be done. 

The cloud discrimination technique performed much better than the earlier tech- 
nique. Evidence from this improvement is the reduction in cloud contamination that 
was contained in the earlier results. This also permited a method of discriminating be- 
tween clouds and aerosols that was not available earlier. 

The sunglint removal improvement added additional data that was previously lost 
to processing. The removal of Rayleigh radiance and the accounting for ozone absorp- 
tion produced a more accurate aerosol particle size index. The estimation of optical 
depth vice albedo gives a better understanding of the radiative transfer characteristics 
of the atmosphere. It also allowed for a correlation between optical depth and APSI 
values which was ver>' useful in identifying the type of aerosol particles. 

The data processed provided several interesting results. A spatial view of the relative 
aerosol size distribution was possible from the APSI values. It was easy to distinguish 
man-made increases over naturally-occuring increases. More evidence was provided that 
showed the ability of the atmosphere to transport aerosol particles great distances, both 
in the North and South Atlantic Ocean as well as in the North and South Pacific Ocean. 

Some evidence was provided for the theory that much of the aerosol particle popu- 
lation, at least in the tropics, is biogenic in origin, and that there may be a biological 
feedback mechanism in the earth's radiation budget. This evidence was provided 
through intercomparing data during an El Nino-Southern Oscillation time period with 
one that came before it. The drastic change, particularly in the eastern North Pacific 
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Ocean, showed that there is some form of climatic influence on aerosol particles in the 
tropics. 

The evidence would have been more substantial if more data could have been proc- 
essed. This may be in the form of continuing to process only April data for a number 
of years and or processing a years worth of data prior to, during and after the El Nino- 
Southern Oscillation event. If the values remain consistent before and after the El 
Nino-Southern Oscillation event and show a decline during the event, this would provide 
more proof of biogenic origins. Additionally, it would be advantageous to have cruise 
data that included aerosol particle collections as well as phytoplankton measurements. 
That, in addition to the satellite data, should be able to prove or disprove this theory'. 

One improvement in the technique that should be implemented in the future would 
be to obtain the APSI from optical depth ratios vice radiance ratios. This should elim- 
inate the slight geometric bias that is observable in Fig. 12. Since the geometric bias 
has been eliminated from the optical depth and the radiance ratio is an approximation 
for the optical depth ratio, it would be better to calculate this directly. This would also 
provide a better correlation between optical depth and APSI values. 

This study presents an excellent tool for the evaluation of aerosol particles on a 
global scale. With a few minor modifications, it hopefully will provide critical data for 
long term global radiation studies. 
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